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The electron localization tensor, Born effective charges, dielectric constants and phonon dispersion 
relations of the skutterudite CoSbs and the filled skutterudite TlFeCo3Sbi2 have been studied using 
density functional perturbation theory. The origin of the low energy peak in the phonon density of 
states of TlFeCo3Sbi2 observed recently by Neutron inelastic scattering [R. P. Hermann et al., Phys. 
Rev. Lett. 90, 135505 (2003)] is attributed to the vibrations of Tl that show only weak coupling 
with the normal modes of the host crystal. Moreover, the dielectric properties of these materials 
show unsual features such as giant Born effective charges and a strong increase of the optical and 
static dielectric tensor in the filled compound. 

PACS numbers: 63.20.Dj,77.22.-d,71.15.Mb 
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I. INTRODUCTION 

CoSba belongs to the class of binary skutteruditesi. 
These materials are characterized by a complex cristalline 
structure containing large voids and four formula units 
per cell. It has been observed that filling the voids with 
electropositive atoms such as Tl or rare earth atoms leads 
to a drastic decrease of the lattice thermal conductiv- 
itji2i^. Combined with good semiconducting properties 
(electrical conductivity, Seebeck coefficient), this makes 
the so-called filled skutterudites potentially interesting 
for thermoelectric applications such as refrigeration or 
energy generation^. 

One mechanism that has been evoked to explain the 
strong reduction of the lattice thermal conductivity in 
filled skutterudites is a strongly anharmonic rattling mo- 
tion of the filling atom. Keppens and co-workers^, used 
inelastic neutron scattering to measure the phonon den- 
sity of states (DOS) of the La atoms in LaFe4Sbi2. They 
found a well defined peak at about 56 cm~^ and a some- 
what broader structure around 121 cm^^ that were as- 
sociated to two La-dominated localized modes. Feldman 
and co-workers^ used a force constant model fitted from 
first-principles calculations to compute the phonon dis- 
persion relations and DOS of LaFe4Sbi2 and CeFe4Sbi2. 
Their results allowed them to refute the hypothesis of an 
anharmonic rare-earth potential and to suggest that the 
two peaks in the La projected DOS are not associated to 
two localized modes but are the results of hybridizations 
between the rare-earth and Sb vibrations. 

In several respects, we can expect the behaviour of Tl 
to be different from that of La and Ce. First, the 4f 
valence states of the rare earth atoms strongly hybridize 
with the transition metal d and Sb p states in the vicinity 
of the Fermi energy^. These hybridizations are absent in 
case of Tl where the 4f states are much lower in energy. 
Second, the electronegativity of Tl is significantly higher 
than that of La and Ce and close to that of Sb. This 
suggests that Tl may have a smaller effect on the electri- 
cal transport properties of CoSbs than the rare-earths. 
Finally, the ionic radius of Tl is significantly larger than 
that of La or Ce. We can therefore expect the dynamical 



properties of Tl to be significantly different from that of 
the rare earths. 

Some evidence of a different behaviour of Tl can be 
found in a recent work of Hermann and co-workersS who 
measured the phonon DOS of various Tl-fiUed skutteru- 
dites. These authors found a single peak arount 40 cm^^, 
in striking difference with the two peaks of the La-filled 
compound. 

Motivated by this workSi, we use first-principles calcu- 
lations to study the lattice dynamics and dielectric prop- 
erties of CoSba and TlFeCo3Sbi2. We show that Tl is 
in a harmonic potential well up to large displacements. 
Moreover, the peak about 40 cm~^ in the phonon DOS 
of the filled compound is due to the vibrational motion 
of Tl that shows only modest hybridizations whith the 
normal modes of the host crystal. Our calculations also 
confirm the existence of a single Tl-dominated peak in 
the DOS of the filled compound. In addition, we observe 
unusual dielectric porperties in CoSba and TlFeCo3Sbi2 
such as giant Born effective charges, or a strong increase 
of the static dielectric constant in the filled compound. 
To our knowledge, this is the first first-principles study 
of Tl-fiUed skutterudites and the first time Born effec- 
tive charges and dielectric constants have been computed 
from first-principles for these compounds. 

Our paper is organized as follows. Sec. ^ contains 
the technical details of our calculations. In Sec. IIIII we 
discuss the structural and electronic properties of CoSba 
and TlFeCo3Sbi2. In Sec. IIVI we discuss their opti- 
cal dielectric constant and effective charges. In Sec. 
we report the zone-center phonons in the filled and un- 
filled compound and deduce their static dielectric tensor 
and infrared refiectivity. Finally, in Sec. IVII we discuss 
the phonon dispersion relations and DOS of CoSbs and 
TlFeCo3Sbi2. 



II. TECHNICAL DETAILS 

Our calculations have been performed within the lo- 
cal density approximation to density functional theory 
(DFT)-^-^-JJ- thanks to the abinit— package. We worked 



at the experimental lattice constant. The atomic posi- 
tions have been relaxed until the forces on the atoms 
were smaller than 5 • lO^'"^ hartree/bohr. For the 
exchange-correlation energy, we used the parametriza- 
tion of Perdew and Wangi^. The valence states were 
computed non-relativistically. The all-electron potentials 
were replaced by norm-conserving pseudopotentials gen- 
erated from scalar relativistic all-electron calculations of 
the free atoms thanks to the FHI98PP codei^. The Sb 
potential has been generated according to the TrouUier- 
Martins schemei^. In order to obtain smooth pseudopo- 
tentials for Fe, Co and Tl we used the Hamann- scheme 
for these elements. The parameters used to generate the 
pseudopotentials are summarized in Tabled The elec- 
tronic wavefunctions have been expanded in a plane wave 
basis. With the pseudopotentials described above a ki- 
netic energy cutoff of 25 hartree was enough to obtain 
well converged results. 



The static ionic charges of Sec. II V Bl have been com- 
puted following Bader's prescription for the partitioning 
of space into atomic bassing-^ as it is described in Ref. 



For the phonon frequencies, a 2 x 2 x 2 grid of special 
k-points^^ was sufficient to obtain well converged results 
while the Born effective charges, dielectric tensors and 
localization tensors required a grid of 4 x 4 x 4 special 
k-points. In fact, for the smaller grid, the violation of 
the charge neutrality for the Born effective charges^^ was 
quite important (0.2 a. u.). Using the larger grid, this 
error could be reduced by a factor of 10 while the phonon 
frequencies changed by less than 3 cm~^. 



The phonon frequencies. Born effective charges, lo- 
calization tensors and dielectric tensors have been com- 
puted from a linear response approach to DFTiSSiSLS^. 
The first-order derivatives of the wavefunctions have been 
calculated by minimizing a variational expression of the 
second-order energy derivatives within the parallel gauge. 
In order to perform the band-by-band decomposition of 
various quantities such as the localization tensor or the 
optical dielectric tensor, these wavefunctions were further 
transformed to the diagonal gauga^^. We computed the 
dynamical matrix explicitely on a 2 x 2 x 2 grid of spe- 
cial q-points. To obtain the full phonon band structure 
and DOS, we used a Fourier interpolation that includes 
the long range dipole-dipole interactions^'*. The lattice 
specific heat has been computed from the DOS as it is 
described in Ref. |2^ This approach is more reliable 
than the one used in the previous studies of Singh and 
co-workeraSi2&. In Refs. |23 and |a, these authors used 
an empirical force constant model that had been adjusted 
throught LDA total energy calculations while in Ref. 7, 
they computed the interatomic force constants for zone- 
center atomic displacements only. 



III. GROUND-STATE PROPERTIES 



A. Crystal structure 

The structure of CoSbs is formed by a bcc lattice [space 
group ImS, lattice constant a = 9.0385 A] with 4 formula 
units per primitive unit cell. It can be described as a dis- 
torted perovskitei with the chemical formula ABO3. In 
CoSbs, the site of the A atoms is empty, the sites of the B 
atoms are occupied by Co and the positions of O by Sb. 
The CoSbg octaedra are tilted so that the Sb atoms form 
rectangular Sb4 rings. The Co and Sb atoms occupy re- 
spectively the 8c (|, J, j) and 24g (0, y, z) Wyckoff posi- 
tions. The relaxed values oiy — 0.33284 and z ~ 0.15965 
arc in exellent agreement with the experimental resultSiS 
of y = 0.33537 and z = 0.15788. 

In the unfilled compound, the 2a (0, 0, 0) Wyckoff po- 
sitions are empty. They are at the centers of large voids 
surrounded by 12 Sb and 8 Co atoms. These voids are 
able to accomodate an electropositive filler such as Tl or 
various rare earth atoms. The electronic configuration of 
Tl is 6s^ 6p*. In order to fill an important fraction x of 
voids, all chemical bonds must be saturated. Therefore, 
one has to compensate the excess charge due to the Tl 
6p electrons by replacing for example one Co per cell by 
Fe or one Sb by Sn^. 

Experimentally, the lattice constant of Tl-fiUcd CoSba 
is found to increase linearly with the amount of filled 
voids'^ for < a: < 0.8. In our calculation, we worked at 
X = 1 and we replaced one Co per cell by Fe. We chose a 
lattice constant of 9.1276 A that corresponds to a linear 
extrapolation of the data of Ref. y to x = 1 . This choice 
is further justified by the fact that the residual stress in 
the filled compound (2.7 • 10^* Ha/bohr'^) is similar to 
the stress in the unfilled compound (2.5 T0~* Ha/bohr'^). 
The relaxed reduced atomic coordinates in TlFeCo3Sbi2 
and CoSba differ by less than 0.003. 

Due to the use of periodic boundary conditions, the Fe 
atoms occupy the same crystallographic site in each cell. 
This situation is somewhat artificial since it neglects any 
disorder introduced by the random substitution of Co 
by Fe in real samples. However, we expect the lack of 
disorder in our calculations to have only modest effects 
on the electronic and dynamical properties of Tl for the 
Fe-compensated compound. It might be more important 
for other compensations such as the substituaion of Sb by 
Sn. First, Fe can only occupy four equivalent sites per cell 
while there are 12 sites that can accomodate Sn. Second, 
the distance between Tl and Sb (3.3897 and 3.3977 A) 
is smaller than the distance between Tl and Co (3.9823 
A). Finally, wc will sec in Sec. IVII that the vibrations of 
Co/Fe are strongly decoupled from the vibrations of Tl 
and located in the high energy region above 200 cm^* of 
the phonon spectrum. At the opposite, the vibrations of 
Sb occupy the region below 200 cm~* and show signifi- 
cant interactions with the vibrations of Tl. 



TABLE I: Parameters used to generate the pseudopotentials of Fe, Co, Sb and Tl. Reported is the type of the pseudopotential 
(H — Hamann, TM = TrouUier-Martins) , the projector that is used as local component, the cutoff radii (bohr) of the s, p, d 
and f channels (in case they have been included in the pseudopotential) and the cutoff radius (bohr) of the partial core density 
(rnic)- The values in brackets indicate the reference energy (eV) that has been used to build the pseudo-wavefunction in case 
we did not use the default value of the FHI98PP code. 



Fe 



Co 



Sb 



Tl 



Type 

Configuration 

Local part 

s-channel 

p-channel 

d-channel 

f-channel 



H 

4s^ 3d^ 

f 

1.2 

1.2 (5.0) 

1.0 

1.6" 

0.9 



4p" 



H 

4si 3d^ 4p" 

f 

1.2 

1.2 (5.0) 

1.0 

1.6" 

0.9 



TM 

5s2 5p^ 

s 

1.8 

2.0 

2.5 (15.0) 

1.5 



6s- 



H 

Sdi" 

P 

1.4 
1.4 
1.4 
1.6" 
1.0 



6p° 



"Default value used by the FHI98PP code 



B. Electronic properties 

The left part of Fig. Q] shows the LDA band structure 
of CoSba between the high symmetry points N, F, and 
P. The top of the valence bands has been fixed at eV. 
There is a good agreement between our results and the 
band structures obtained in previous studiesSSiS^iSMi. 
We observe the presence of two well separated groups 
of valence bands. The bands between -13 eV and -7 eV 
are mainly composed of Sb 5s states while the valence 
bands above -6 eV result from hybridizations between 
Sb 5p and Co 3d bands. 

There is some controversy about the energy of the band 
gap. Eg, in CoSba. Experimentally, values between 0.03 
eV and 0.7 eV have been reported [see Refs. .32^331341 
and references therein]. From a theoretical point of view, 
there is a similar strong discrepancy between the values 
in the litterature. The bandgap is defined by a highly dis- 
persive band at the F-point. We obtain a value of 0.24 eV 
in reasonable agreement with other first-principles LDA 
and GGA calculations [0.05 eV^Sa, 0.22 eViS^, 0.195 & 
0.330 cV^", 0.140 eViSi]. This discrepancy between the 
theoretical results has been attributed to a strong depen- 
dence of Eg on the position of the Sb atoms together with 
the well-known problem of the LDA and GGA to predict 
the correct bandgap of semiconductor a^^i^" . 

In the TlFeCo3Sbi2, [right part of Fig. Q], the bandgap 
(0.20 eV) is slightly smaller than in the unfilled com- 
pound. In addition, we observe two additional groups of 
bands that are related to the atomic orbitals of Tl. The 
Tl 5d states form a flat set of bands that are in the same 
energy region as the Sb 5s states. Due to the lack of 
dispersion, these states can be considered as inerte. At 
the opposite, the Tl 6s states form a band about -7 eV 
that shows a significant dispersion related to a covalent 
interaction with the electronic states of the host crystal. 
This point will be investigated in more detail below. 

In some cases, the filling of the voids in skutterudites 
does not only affect the lattice thermal conductivity but 
also the electronic transport properties. For example, the 



CoSb^ 



TIFeCogSbig 




FIG. 1: Electronic band structure in CoSbs and 

TlFeCogSbia. 



Ce 4f states in CeFe4Sbi2 hybridize with the Fe 3d and 
Sb 5p states to form a set of flat conduction bands in 
the vincinity of the Fermi level^. Due to the low disper- 
sion of these bands, the electron effective mass in this 
compound is quite high implying a low electron mobil- 



ity 
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As it can be seen in the right part of Fig. ^ 



the effect of Tl filling on the electronic band structure of 
CoSbs is less important. The only additional bands ly 
well below the Fermi level and will therefore not affect 
the electronic transport properties. However, in order to 
conclude definitively about this point, it might be impor- 
tant to perform a more intensive study of the electronic 
band structure of TlFeCo3Sbi2 that eventually takes into 
account spin-orbit coupling. 

In order to get additional informations on the elec- 
tronic properties of CoSbs and TlFeCo3Sbi2, we com- 
puted the localization tensor (see Ref. ,2fl and references 
therein) and the electronic density for selected states 
in these compounds. As we have shown in a previous 



work^^, the decomposition of the locahzation tensor into 
contributions of individual groups of bands can act as 
a sensitive probe to study hybridizations within a sohd. 
For example, a state that has the same variance in a solid 
and in the isolated atom can be considered as chemically 
inerte. At the opposite, a modification of the variance is 
usually attributed to covalent interactions of the corre- 
sponding orbitals. 

In CoSbs and TlFeCo3Sbi2 the localization tensor is 
isotropic. In Table Ull we report the variances of the two 
(resp. three) groups of bands defined in Sec. IIIIBI as well 
as the variances of the atomic orbitals Sb 5s, Tl 5d and Tl 
6s. In case of Tl, we used linear response calculations to 
obtain the variance of the 6s state of an isolated T1+ ion. 
We adopted this configuration because it is close to the 
configuration we expect for Tl in TlFeCo3Sbi2. More- 
over, it allowed us to avoid some problems related to the 
linear response calculation of the localization tensor for 
partially occupied states. To use the same approach for 
Sb, we had to chose Sb^+ as reference configuration. For 
the 5s state, we obtained a value of 1.51 bohr^. Unfor- 
tunately, the character of the bonding of Sb in CoSba is 
rather covalent than ionic. The Sb'^"'" configuration seems 
therefore somewhat artificial. To compute the variance 
for the neutral Sb atom, we used the FIII98PP— code to 
compute the all-electron wavefunction of the 5s state and 
its variance as a matrix element of the squared position 
operator x^ . 



in case the covariance between Sb 5s and Tl 5d bands is 
not zero since this quantity is necessarily negative. 

A stronger variation is found for the variance of 
the highest occupied bands (Sb 5p -I- Co 3d). In 
TlFeCo3Sbi2, these electrons are more delocallized than 
in CoSbg. 

The most spectacular evolution undergoes the variance 
of the Tl 6s state. In the solid, this quantity is about 9 
times larger than in the isolated ion. This result suggests 
a significant interaction between the Tl 6s orbitals and 
the electronic states of the host crystal. To get more in- 
formations about these hybridizations we show in Figure 
13 the electronic pseudo-density associated to the Tl 6s 
band. As can be seen, these electrons are partially delo- 
callized on the Sb atoms of the cage. This result suggests 
that there is a significant covalent interaction between Tl 
and Sb due to the Tl 6s electrons. Such an interaction 
would be coherent with the fact that Tl can be inserted 
into the voids of the host crystal in spite of the fact that 
its ionic radius (1.84 A) is larger than that of the rare- 
earth fillers Ce (1.28 A) and La (1.50 A) and that it is 
close to the radius of the cage (about 1.892 A2I). 

In order to obtain further insight on the electronic 
properties of skutterudites and the effect of different fill- 
ing atoms on the band structure it might be interesting to 
build generalized Wannier functions^S*^ for these com- 
pounds. 



TABLE II: Variances {r'^)c (Bohr^) of the groups of bands 
in CoSbs and TIFeCo3Sbi2 and of the Sb 5s, Tl 5d and Tl 6s 
electrons calculated on an isolated Sb atom and Tl^ ion. 



Bands 


CoSbg 


TlFeCogSbia 


Atom 


Sb 5s (+ Tl 5d) 

T16s 

Sb 5p + Co 3d 


3.06 
5.11 


2.63 
18.06 
5.85 


1.86", 0.85" 
2.14 


E, (eV) 


0.24 


0.20 


— 



"Sb 5s orbital of an isolated Sb atom. 
"Tl 5d orbitals of an isolated T1+ ion. 



In CoSba, the variance of the Sb 5s bands is larger 
than the variance of the corresponding atomic orbital. 
This result confirms the participation of these electrons 
in the Sb-Sb bonding observed previously—. In the filled 
compound, the Tl 5d states occupy the same energy re- 
gion as the Sb 5s states. Therefore, the formalism of 
ReL [22 cannot be applied to study the localization of 
these bands individually. Nevertheless, we can obtain 
some information on the modification of the Sb 5s states 
in TlFeCo3Sbi2. Due to the flat dispersion of the Tl 5d 
bands observed in Fig. ^ it is reasonable to assume that 
their variance in the solid is close to the value in the iso- 
lated T1+ ion. If we neglect the covariance between Sb 5s 
and Tl 5d bands, we obtain a value of 3.37 Bohr^ for the 
variance of the Sb 5s bands in TlFeCo3Sbi2. This result 
suggests that the Sb 5s electrons are more delocallized in 
the filled compound than in the unfilled. It remains true 




FIG. 2: Electronic pseudo-density of the Tl 6s band in 
TlFeCo3Sbi2. The isodensity value corresponds to 10 % of 
the maximum density of this state (0.0013 electron/bohr'^). 
Tl (not shown) occupies the center of the cage. The black 
and grey spheres represent respectively Sb and Co/Fe. This 
figure has been realized thanks to the XCrySDen— crystalline 
and molecular structure visualisation program. 

In summary, the main differences between the elec- 



tronic band structures of CoSbs and TlFeCo3Sbi2 affect 
the region well below the Fermi level. On the one hand, 
the Tl 5d states show a flat dispersion suggesting that 
these states are chemically inerte. On the other hand, 
the Tl 6s states show a significant dispersion related to 
a partial delocalization of the corresponding electrons on 
the Sb atoms of the cage. These hybridizations are ex- 
pected to affect the dielectric and dynamic properties of 
Tl in the filled compound. 



IV. DIELECTRIC PROPERTIES 
A. Optical dielectric tensor 

In CoSba, the dielectric tensor is isotropic. The theo- 
retical value £00 = 31.67 overestimates the experimental 
resull^i of 25.6 by about 20 % as it is typical in LDA 
calculations^. Similar large dielectric constants as in 
CoSba have been computed for other narrow bandgap 
semiconductors such as AI2RU [principal values: 18.9, 
22.9, 20.7)^. Moreover, they have been measured in a 
series of skutterudites such as C0AS3 (26.25)^, UFe4Pi2 
(17) or CeFe4Pi2 (31)^. 

In the filled compound, the dielectric tensor is no more 
isotropic but has weak off-diagonal elements. These ele- 
ments are about four orders of magnitude smaller than 
the diagonal ones and we will not consider them in the 
discussion that follows. We obtain a value of 38.57 for the 
dielectric constant in TlFeCo3Sbi2 that is significantly 
larger than in unfilled CoSba . Unfortunatley, we are not 
aware of any experimental data to which we can com- 
pare this result. In spite of the error on its absolute 
value, we expect the evolution of the dielectric constant 
to be correctly reproduced by our calculation since it is 
reasonable to assume that the LDA error on Eoo is of the 
same magnitude in both compounds. 

The increase of £00 can have distinct origins. On the 
one hand, based on a simple Clausius-Mossotti model—, 
the polarizabilities of the T1+ ions add to the dielectric 
susceptibility of the unfilled compound to yield a larger 
£00 in TlFeCo3Sbi2. On the other hand, the increase can 
be due to the dielectric constant of the host crystal that 
is increased due to the presence of Tl and Fe in the filled 
compound. To check which mechanism applies in CoSbs, 
we decomposed the optical dielectric tensor into contri- 
butions originating from the individual groups of bands 
of Fig. n This decomposition reveals that the deeper 
bands (Sb 5s, Tl 5d and Tl 6s) are only weakly polar- 
izable by an electric field and contribute less than than 
0.1 % to £00 and that the contribution of the Sb 5p and 
transition metal 3d electrons is strongly increased in the 
filled compound. We therefore propose that the observed 
increase of £00 is due to a modification of the dielectric 
constant of the host crystal while the polarizability of 
Tl''" plays only a minor role. 



B. Effective ionic charges 



In Table IIIII we report the static ionic charges Qi 
in CoSba and TlFeCo3Sbi2 obtained from a topological 
analysis following Bader's prescription for the partition- 
ing of space into atomic bassins— . In this approach, one 
has to determine surfaces that obey the zero-flux condi- 
tion for the electron density n(r): Vn ■ N — where TV 
is the vector normal to the surfaces. The integration of 
n(r) within the atomic bassins yields the ionic charge Qj 
of each topological atom i. 

In both compounds the chemical bonds have a negli- 
gible ionic character as it is revealed by the small Q^. 
At the opposite, to these small static charges, the dy- 
namical Born effective charges Z* (Table IIV|I are very 
large, comparable to the giant effective charges in ferro- 
electric ABO3 compounds— and other narrow bandgap 
semiconductors such as Al2Ru^. In the ferroelectrics, 
the amplitude of the Z* can be explained thanks to the 
Harrison model^. A relative displacement of the B and 
O atoms generates a giant dipole moment due to a dy- 
namical change of hybridization between O 2p and B d 
atomic orbitals^. In CoSba, the highest valence bands 
are basically due to hybridizations between Co 3d and 
Sb 5p atomic orbitals^. Moreover, the structure of this 
compound is formed of distorted CoSbg octahedra simi- 
lar to the BOg octahedra of the ABO3 ferroelectrics. It 
seems therefore plausible to assume that a similar mech- 
nisme takes place in CoSb3 and that the large effective 
charges are generated by dynamic changes of hybridiza- 
tions between Co 3d and Sb 5p atomic orbitals. 

In skutterudites, giant effective charges have been mea- 
sured from infrared spectroscopic for UFe4Pi2 {\ZpJ = 
8.9) and CeFe4Pi2 {\Z*pJ = 11)^. In CoSb3, we obtain^ 
a value of 6.04 for \Z'^J from the amplitude of the LO- 
TO splitting reported in Ref. 51 and the optical dielectric 
constant of Ref. |4l| This value is in good agreement with 
the results reported in Tab. IIVI 



TABLE III: 
TlFeCogSbia, 



Bader ionic charges Qi (|e |) in CoSba and 





CoSba 




TlFcCosSbia 






Tl 


0.29 


Coi 


-0.53 


Fe 


-0.36 


C02 


-0.53 


C02 


-0.48 


Sb 


0.17 


Sb 


0.12 



The effective charge of the filling Tl atom is also highly 
anomalous. It is significantly larger than its static ionic 
charge reported in Table UTTl Moreover, it is larger than 
its nominal ionic charge of -1-1 that would be expected 
if there was a complete transfer of the 6p electron to 
the atoms of the cage. Giant effective charges have been 
reported previouslj*^ for the filling atoms U and Ce in 



FeP4. Their values of Z, 



u 



8.9 and Z, 



Ce 



8.7 are even 
larger than the value of Z^^ in CoSb3. 

The effective charges of Co and Sb are similar in the 



TABLE IV: 


Born 


effective 


cliarges 


(|e"l) 


in CoSba and 


TlFeCosSbia 


















CoSba 






TiFeCoaSb 


L2 










Tl 


3.64 
-0.02 
-0.05 


-0.05 

3.64 

-0.02 


-0.02 

-0.05 

3.64 


Coi 


-6.89 


0.82 


0.34 


Fe 


-9.83 


1.13 


-0.20 




0.34 


-6.89 


0.82 




-0.20 


-9.83 


1.13 




0.82 


0.34 


-6.89 




1.13 


-0.20 


-9.83 


C02 


-6.89 


-0.82 


-0.34 


C02 


-7.28 


-0.75 


0.08 




-0.34 


-6.89 


0.82 




-0.09 


-7.30 


0.73 




-0.82 


0.34 


-6.89 




-0.64 


-0.10 


-7.34 


Sb 


2.63 


0.00 


0.42 


Sb 


2.43 


-0.10 


-0.04 




0.00 


1.94 


0.00 




-0.03 


1.87 


0.01 




0.68 


0.00 


2.32 




1.01 


0.12 


2.13 



filled and the unfilled compound. The excess effective 
charge introduced by Tl can therefore not be compen- 
sated by these two atoms in order to satisfy the charge 
neutrality conditionSS ^^ Z* ^^ — 0. As can be seen in 
Table HVI \Z*p\^ is about 3|e~| larger than the absolute 
value of the effective charge of the Co atom that it substi- 



tutes. As mentioned in Sec. IIII Al the 6p electron of Tl 
is compensated by a hole in the 3d orbitals of Fe. In case 
of the Born effective charges, Fe plays a similar role: due 
to the fact that \Z*p^ is larger than |.^5oIj the anoma- 
lous effective charge of Tl is mainly compensated by the 
anomalous effective charge of Fe. A plausible mechanism 
for this compensation could be that the Tl 6p electron 
is partially shared between Tl and Fe so that it might 
be responsible for a transfer of charge during a relative 
displacement of these two atoms. 



been suggested previously"^^. Close to the origin, the 
energy changes quadratically with the atomic positons 
[dashed line] . The whole energy region of Fig. O can 
be fit accurately by a polynome of degree six [full line]. 
The harmonic frequency is 52 cm~^. This so-called "bare 
frequency" is quite independent of the direction of the 
atomic displacement. For the directions [1 1 0] and [1 1 
1], we obtaiend values that differ by less than 0.3 cm~^ 
from the one of Fig. In case of Tl, the bare fre- 
quency is smaller than the one of Ce (68 cm~^) and La 
(74 cm^^) in FeSb4^. This difference is attributed to 
the larger mass of Tl and to the curvature of the po- 
tential energy at the origin that is significantly smaller 
in case of Tl (2.0371 eV/A) than for Ce (2.3943 eV/A) 
and La (2.8005 eV/A). This result is surprising. On the 
basis of the crystal radii^^ ^f rpii+ {^^M A), Ce4+ (1.28 
A) and La'^^ (1.50 A) in a 12-coordinate site we could 
expect the repulsion bctwenn ionic cores to be more im- 
portant between Tl and Sb than between Ce/La and Sb. 
The eneregy should therefore increase more rapidly for 
an off-center displacement of Tl and the curvature of the 
potential energy well should be stronger. However, in 
Sec. IIII Bl we saw that there is a significant covalent in- 
teraction between Tl and Sb. The simple argument of 
the crystal radii can therefore not be applied in case of 
TlFeCo3Sbi2 because the atoms do not behave like rigid 
spheres. 

The large ADP's of Tl have been attributed to an 
Einstein oscillator— with a characteristic frequency of 36 
cm~^. This value is close to the Einstein frequencies 
deduced from specific heat measurements and inelastic 
neutron scattering experiments^ (39 ± 1 cm~^) but it is 
significantly lower than the computed Tl bare frequency. 
This difference between the bare frequency and the nor- 
mal mode frequency observed experimentally will be fur- 
ther discussed in Sec. IV Bl 



V. ZONE-CENTER PHONONS 

A. Potential energy 

To understand the effect of Tl-fiUing on the lattice 
thermal conductivity of CoSba, it is mandatory to in- 
vestigate the dynamics of these atoms in their cages. 
Experimentally, unusually large atomic displacement pa- 
rameters (ADP's) have been observed for various filling 
atom o?'^^ . To decide whether these ADP's are due to 
quasi-harmonic vibrations about the center of the cage 
or to the random hopping between off-center siteE^=^», it is 
important to investigate the shape of the potential energy 
well. 

To determine the shape of the Tl potential, we com- 
puted the variations of the total energy for atomic dis- 
placements along various directions. Figure 13 shows the 
results for the cristallographic direction [10 0]. We see 
that the center of the cage is the global minimum of the 
potential energy well and not a saddle point as it has 
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FIG. 3: Potential energy associated to a cooperative displace- 
ment of Tl along [10 0]. 



B. Normal mode frequencies 

The unfilled skutterudite CoSbs belongs to the space 
group Im3. At the F-point, the zone center optical 
phonons can be classified according to its irreductible 
representations into 



2A„ + 2E„ + AFq + 2Au + 2E^ + 7F„. 



(1) 



Only the F„-modes are infrared active. At the F-point, 
they are split into transverse (TO) and longitudinal (LO) 
modes. The Ag, Eg and Fg are Raman active. 

In Table Ivl we report the frequencies (cm^^) of the 
TO and LO modes in CoSbg and TlFeCogSbia. The val- 
ues in the unfilled compound are compared to the first- 
principles and force-constant (PC) model calculations of 
Feldman and Singh [model A of Ref. ,26] and to the 
experimental data obtained by infrared^ and RamanS 
spectroscopy. The overall agreement between theory and 
experiment is quite good. Especially the phonon frequen- 
cies in the low energy region — the region of the Tl vibra- 
tions in the filled compound — are accurately reproduced 
by our first-principles calculations. 

We should note that the Raman measurements of Ref. 
|53 were performed on polycrystalline samples. As a con- 
sequence, only the Ag modes were identified in this study 
while the remaining Raman lines can be either associated 
to Eg or Fg modes. In TablelVI we identify four of the ad- 
ditional frequencies to the theoretical Fg modes. In fact, 
the experimental frequencies of 82, 109 and 152 cm~^ are 
close to the three lowest Fg modes obtained in our study 
and by Feldman and Singh^^. In case of the mode at 178 
cm~^, this identification is less clear since its frequency 
is close to the highest Eg mode and the highest Fg mode. 
In addition to these four modes, the authors of Ref. 53 
report a mode at 60 cm^^ [not reported in Table Iv] that 
is not confirmed by our calculation neither by the cal- 
culation of Feldman and Singh ^^. Therefore, we do not 
associate this Raman line to a phonon mode in CoSba. 

In spite of the giant effective charges, the transverse 
and longitudinal F„ modes are rather close in energy. 
This small amplitude of the LO-TO splitting can be at- 
tributed to the large optical dielectric constant (see Sec. 
IIV A|) that effectively screens the macroscopic electric 
field generated by the longitudinal polar lattice vibra- 
tions. At lower energies, this screening is almost perfect 
and the LO and TO modes are nearly degenerated. At 
higher energies, we observe a small increase in the fre- 
quency of the LO modes with respect to the TO modes. 
We will see in the next section that the low and high 
energy regions are respectively dominated by vibrations 
of Sb and Co. The larger amplitude of the LO-TO split- 
ting at higher energies can therefore be attributed to the 
giant effective charges of Co that are about three times 
larger than those of Sb. 

In the filled compound TlFeCo3Sbi2, the presence of 
Fe breaks the symmetry of the crystal lattice. As a con- 
sequence, it is no more possible to classify the phonon 



modes according to Eq. (^. Nevertheless, the eigenvec- 
tors of most phonons are only slightly modified in the 
filled compound. By computing their overlap, it is there- 
fore possible to associate each phonon in TlFeCo3Sbi2 to 
a mode in the unfilled compound. Due to the symme- 
try breaking, the degeneracy of some phonons is left in 
TlFeCo3Sbi2. In this case, we report the frequencies of 
all the resulting modes if they differ by more than one 
cm"^. The first and second values in the last column of 
Table IVI are respectively one and two times degenerated. 

The TO modes that have the strongest overlap with 
the i^^f modes of the unfilled compound, are polar. The 
single and double degenerate modes are respectively po- 
larized in the [1 1 1] direction and in the plane perpen- 
dicular to [1 1 1]. The frequencies of the corresponding 
LO modes are reported at the end of Tab. In addition 
to the Fu modes other modes become polar in the filled 
compound. Among them, the highest Au and E^ modes 
acquire the strongest polarity. Nevertheless, the splitting 
between the corresponding TO and LO modes is rather 
weak and not reported in Table IVI 

At 40 cm~^, we observe three additional polar modes 
that are absent in the unfilled compound. Experi- 
mentally, a peak in the phonon density of states of 
TlFeCo3Sbi2 has been observed at this frequency^. This 
peak has been attributed to a locallized vibration of Tl in 
his cage. Our calculation confirms the strong contribu- 
tion of Tl to these modes. About 88 % of the eigenvectors 
are due to the vibrations of Tl. The frequencies of these 
modes are significantly smaller than the bare frequency 
of a pure (100 %) Tl vibration reported in Sec. IV Al 
(52 cm~^). This difference between the bare frequency 
and the frequency of the local Tl mode (LM) can be at- 
tributed to hybridizations between the pure Tl vibrations 
and the phonons of the unfilled compound. 

To the authors knowledge, no experimental data on 
the zone-center phonons of TlFeCo3Sbi2 are available. 
However, in a recent work, Feldman and co-workers^ 
measured the frequencies of the Raman active modes in 
Lao.75Fe3CoSbi2. In addition, they used first-principles 
calculations to compute the frequencies of all zone-center 
optical phonons in LaFe4Sbi2. The results of teir studies 
are reported in Table IVI where they are compared to our 
theoretical results for TlFeCo3Sbi2 although they have 
been obtained for a different compound. However, the 
Raman active modes are dominated by the dynamics of 
Sb. In Sec. IVII we will see that the vibrations of Sb 
are almost decoupled of Co/Fe. We can therefore expect 
the frequencies of the Raman active modes to be similar 
in TlFeCo3Sbi2 and the La-filled skutterudites. In Ta- 
ble we see that this is indeed the case and that the 
agreement between our results and those of Feldman and 
co-workers is quite good for most of the Sb-dominated 
modes between 80 and 200 cm~^. The frequency of the 
La-dominated mode in LaFe4Sbi2 is larger than that of 
the Tl-dominated mode in TlFeCo3Sbi2. This result will 
be discussed more in detail in Sec. IVII 



TABLE V: Frequencies (cm~^) of the transverse and longitu- 
dinal zone-center optical phonons in CoSbs and TlFeCo3Sbi2. 
The phonons in the unfilled compound are classified according 
to the irreductible representations of the space groupe ImS. 
Each phonon in TlFeCo3Sbi2 is associated to the mode of the 
unfilled compound for which the overlap between eigenvectors 
is the strongest. In case the degeneracy of a mode is left in 
the filled compound, we report the frequencies of all resulting 
phonons. The first and second values in column 6 are respec- 
tively one and two times degenerated. The values of Ref. y in 
the two last columns have been obtained for La-filled FeSbs 
for which the lattice dynamics is expected to be similar to 
that in CoSbs. 
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C. Infrared reflectivity and static dielectric tensor 

The static dielectric constant Eq can be decomposed 
into an electronic and an ionic par^ 



EO — £0 






Ae. 



(2) 



Asm represents the contribution of an individual zone- 
center TO mode and can be computed from the infrared 
oscillator strenght Sm-, the phonon frequency ujm and the 
unit cell volume Qq 



A£„,, = 



An S„ 



(3) 



In Table IVll we report the decomposition of Eq in 
CoSbs and TlFeCo3Sbi2. In the imfilled compound, 
the ionic contribution is dominated by the F„ mode 
at 232 cm^^ that has the strongest oscillator strength 
{Sm = 14.29 • lO^^a.u.). Such a large value can be com- 
pared to the oscillator strenghts of the high polar modes 
in ferroelectric oxides such as LiNbOgS. In addition, 
the frequency of this mode is of the same magnitude 
as the frequency of the most polar modes in LiNbOs. 
However, the Ae^ is small compared to the ionic con- 
tribution to eo in the ferroelectrics. This surprising re- 
sult is due to the fact that CoSba has a relatively open 
structure characterized by a large unit cell volume. In 
CoSbs, the primitive unit cell contains 16 atoms and has 
a volume fio = 369. 20A'^ that is more than three times 
larger than the volume of the 10 atom unit cell in LiNbOa 
{^0 = 101.70A3). 



TABLE VI: Electronic and ionic contributions of individual 
phonon modes to the static dielectric constant in CoSba and 
TlFeCoaSbia. 



CoSba 

Lu,„ (cm~-^) Ae,, 



TlFeCo3Sbi2 

LJm (cm~-^) ASn 



Electronic 



31.67 



38.57 
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Au 
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F„ 
F„ 
F„ 
Au 
Eu 
F„ 
F„ 
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0.02 
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0.38 
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0.23 
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0.09 
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0.13 


232 


6.44 


217,215 


8.04 






224 


0.93 






228 


0.06 


245 


0.49 


248,233 


0.94 


258 


1.52 


238,261 


2.07 



Total 



42.44 



63.83 



In TlFcCo3Sbi2, the static dielectric tensor is no more 
isotropic but presents weak off-diagonal elements. As 
for the optical dielectric tensor, these elements are about 
two orders of magnitude smaller than the diagonal ones. 
To a good approximation, we can therefore neglect them 
and characterize the dielectric tensor by its single diag- 
onal element eo- The decomposition of eo in the filled 
compound is more complicated than in the unfilled. The 
values reported for the F„ modes in Tab. IVII correspond 
to the sum of the Ae„i of the three modes that have been 
splitted in the filled compound (see Sec. IV B|) . 

In TlFeCo3Sbi2, eo is about 50 % larger than in un- 
filled CoSbs. Part of this evolution can be attributed 



to the electronic contribution that increases by 6.89 as 
discussed in Sec. IIV Al and to the ionic contribution of 
the F„ modes that is shghtly larger in the filled com- 
pound. In addition, due to its lower symmetry, the A„ 
and F„ modes become infrared active in the filled com- 
pound although thcyr contribution to eo is quite weak. 
The most important part (10.77) comes from the Tl 
dominated mode at 42 cm~^. It's oscillator strength 
{Sm — 0.79 • 10~'^a.u.) is rather low compared to that 
of the Co/Fe dominated modes at high energy. But due 
to its low frequency, the constribution of this additional 
mode in the filled compound is very important. This 
result demonstrates that Tl-filling has an important in- 
fluence on the static dielectric constant of CoSbs. It sug- 
gests that materials with a structure characterized by 
large empty voids might be interesting candidates for ap- 
plications that require a well defined dielectric constant 
since the value of eo can be tuned by filling the voids (or 
a fraction of voids) with foreign atoms. 

From the analysis of the results presented in this paper, 
we can identify some general criteria that should be satis- 
fied by an element in order to be efficient for such a tuning 
of the dielectric constant. First, the Born effective charge 
of the filling atom should be as high as possible. This re- 
quires a certain amount of hybridizations between the 
foreign atoms and the atoms of the host crystal. Second, 
the bare frequency of the filling atom should be as low as 
possible since we expect the normal mode frequency to 
increase with the bare frequency. On the one hand, this 
criterion requires that, in spite of the hybridizations nec- 
essary for a high polarizability and Born effective charge, 
the filling atom should only be weakly bound to ensure a 
low curvature of the potential energy well. On the other 
hand, the mass of the filling atom M should be high since 
the bare frequency is proportional to 1/yM. 

Fig. 2] shows the infrared reflectivity of CoSba (a) 
and TlFeCo3Sbi2 (b) computed from the Born effective 
charges and phonon frequencies and eigenvectors as it is 
described in Ref. .54. Since this approach neglects the 
damping of the phonon modes, the curves saturate to 1. 
For the unfilled compound, there is a qualitative good 
agreement between Fig. 0| (a) and the reflectivity mea- 
sured by Lutz and Kliche^i*^. The most striking differ- 
ence between the reflectivity of CoSba and TlFeCo3Sbi2 
is the presence of a well defined band below 50 cm"^ 
that is related to the Tl-dominated mode and the large 

In addition, we ob- 



Tl Born effective charge (Sec. IIVB|) 
serve the appearance of a few lines associated to the A^ 
and E„ modes that become infrared active in the filled 
compound as discussed above. Unfortunately, we are 
not aware of any experimental refiectivity measurements 
for TlFeCo3Sbi2. However, Dordevic and co-workers'*^ 
measured the refiectivity of MFe4Pi2 filled skutterudites 
(M — La, Th, Ce, U). They observed a similar strong 
band associated to the M-dominated mode in these com- 
pounds. This shows that infrared spectroscopy is an ap- 
propriate technique to study the dynamics of the filling 
atoms in skutterudites. 




FIG. 4: 

(b). 



Frequency (cm" ) 



Infrared reflectivity of CoSba (a) and TlFeCo3Sbi2 



VI. PHONON DISPERSION CURVES, DENSITY 
OF STATES AND SPECIFIC HEAT 



Fig. 13 shows the phonon dispersion relation and pro- 
jected DOS of CoSba (a) and TlFeCo3Sbi2 (b). In both 
compounds, the vibrations of Co/Fe and Sb are mainly 
decoupled and occupy respectively the high energy re- 
gion above 220 cm~^ and the low energy region below 
200 cm-i. 

The computed DOS of CoSbs is in good agreement 
with the DOS obtained from various force constant model 
calculations by Feldman and SinghSi.. At low frequency, 
the agreement is nearly perfect. All calculations also pre- 
dict a large gap at about 200 cm~^ and a smaller gap at 
about 100 cm~^. We do not properly reproduce the two 
well defined peaks above 200 cm~* reported in Ref. l26l 
(calculations) and Ref. |5a (experiment). This might be 
due to a lack of accuracy in our calculation of the high 
energy part of the spectrum. In Tab. we see that the 
frequencies above 200 cm~^ tend to underestimate the 
experiment while those below 200 cm~^ are more accu- 
rate. The small deviations at high energy do not affect 
the following discussions that focus on the lower part of 
the spectrum. 

In the filled compound, an additional branch at about 
40 cm~* gives rise to a well defined peak in the DOS that 
is not present in the unfilled compound. A peak at the 
same frequency has been observed by inelastic neutron 
scattering^ in the DOS of Tl-fiUed CoSba. By comparing 
the DOS of the filled and unfilled compound, the authors 
of Ref. 9 attributed this peak to a locallized oscillation 
of Tl. As can be seen in Fig. [3 the peak at 40 cm^* 
is strongly dominated by the vibrations of Tl that show 
only weak hybridizations with Sb. This result gives there- 
fore a strong argument in favour of the interpretation of 
Ref. la 
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Experimentally, a difference has been observed be- 
tween the DOS of Tl- and La-filled skutterudites. The 
authors of Ref. |3 report the existence of a single Tl- 
dominated peak in TlFeCo3Sbi2 at 40 cni~^ while the 
authos of Ref. 5 claim that there are two La-dominated 
peaks in the DOS of LaFe4Sbi2 at 56 and 121 cm^^. 
Our calculation clearly confirms the existence of a sin- 
gle peak in TlFeCo3Sbi2 (sofid line in Fig. EJ. In Sec. 
IV Al we saw that the bare frequency of Tl (52 cm~^) is 
significantly smaller than that of La (74 cm~^), a differ- 
ence that must be attributed to the larger Tl mass and 
a weaker curvature of the Tl potential energy well. If we 
redo the calculations for a fictitious compound in which 
we artificially reduce the mass of Tl to that of La, the Tl 
vibrational frequency is shifted to 48 cm~^ and yields a 
marginal increase in the DOS around 100 cm~^ (dashed 
line in Fig. 0). Further, if we artificially strengthen the 
interaction of Tl with the lattice by increasing the inter- 
atomic force constants in order to fit the curvature of the 
potential energy well of La, the main Tl vibrational peak 
is shifted at 51 cm~^ whereas we observe the appearance 
of a second peak around 100 cm^^. This strongly sug- 
gests that the peak at 121 cm"-'^ reported by Keppens and 
co-workers5. cannot be attributed to an additional local- 
ized mode but rather to phonon modes that are formed 
by coupled Sb-La vibrations. Generally, it appears that 
the coupling between the La and Sb vibrations is more 
important than between those of Tl and Sb because in 
the former case, the bare frequency is closer to the en- 
ergy region of the Sb-dominated optical phonons. This 
interpretation of the second peak is consistent with the 
one drawn by Feldman and co-workers in Ref. |a- 

As discussed above, the vibrations of Co/Fe are lo- 
cated in the high energy region of the phonon spectrum 
so that they are mainly decoupled from the vibrations 
of Tl. This result allows us to justify the approximation 
used in our study to negelect any disorder due to the ran- 
dom substitution of one Co atom per cell by Fe (see Sec. 
HTj) . Since the vibrations of Tl and Co/Fe are only weakly 
coupled, we expect the infiuence of the disorder on the 
lattice dynamics of Tl to be small. This hypothesis is 
further confirmed by the observation that the Tl-peak 
in the experimental DOS^ is broader in the Sn compen- 
sated compound TlCo4SbiiSn than in TlFeCo3Sbi2. In 
fact, the vibrations of Tl are located in the same en- 
ergy region as that of Sb/Sn. As a consequence, there 
are significant hybridizations between the vibrations of 
those atoms. We expect therefore that the disorder due 
to the random substitution of one Sb per cell by Sn has 
a stronger influence on the dynamics of the filling atom 
than that of Co by Fe. 

Due to the rather flat dispersion of the optical 
branches, we expect the lattice tcrmal conductivity of 
CoSba to be dominated by the acoustic phonons. As 
can be seen in Fig. [S] (b), the Tl-derived modes cut 
through the acoustical branches. A similar effect has 
been observed recently in Sr-fiUed Ge clathrates from 
atomic simulations making use of empirical interatomic 



potentials^. Using molecular dynamics simulations, the 
authors of Ref. |53 showed that the Sr derived modes 
reduce the lattice thermal conductivity of the clathrates 
by approximately one order of magnitude by scattering 
the heat-carrying acoustical phonons. This scattering is 
enhanced at the frequency of the Sr- modes because of the 
resonant interaction with the acoustical phonons^SiSiSS. 
It is therefore reasonable to assume that the resonant 
interaction between the Tl-derived modes and the acous- 
tical branches has a similar strong effect on the thermal 
conductivity of CoSbs and that it is one major source of 
the important decrease observed experimentally^. 

Fig. [3 shows the heat capacity at constant volume 
of CoSba and TlFeCo3Sbi2 computed by integrating the 
phonon DOS as described in Ref. ■25'. Our results are 
compared to the experimental values of Ref. j^ measured 
for CoSba and Tlo.8Co4SbiiSn at constant pressure. For 
the unfilled compound, the agreement between theory 
and experiment is excellent. In the filled compound, we 
observe a slight deviation from the experimental data. 
Part of the discrepancy can probably be attributed to 
the fact that the experiment has been performed on a Sn- 
compensated sample with only partial void-filling. Espe- 
cially the substitution of Sb by Sn is likely to have a 
stronger effect on the specific heat at low temperatures 
than the substitution of Co by Fe since the Sb/Sn vibra- 
tions are at lower energy than that of Co/Fe. 



VII. CONCLUSIONS 

In this work, we used first-principles DFT calculations 
to study the electronic, dielectric and dynamical proper- 
ties of CoSba and TlFeCo3Sbi2. 

Because of hybridizations between Tl and Sb, the Tl 6s 
electrons are partially delocalized on the Sb atoms. This 
result is consistent with the fact that Tl can be inserted 
into the voids of CoSba inspite of its rather large ionic 
radius. 

The electronic dielectric constant of CoSba and 
TlFeCo3Sbi2 are quite large with a significantly larger 
value in the filled compound. This increase cannot be ex- 
plained in the framework of a Clausius-Mossotti model 
where the T1+ polarizability is supposed to add to the 
dielectric constant of CoSba. Instead, it is the dielectric 
constant of the host crystal itself that is increased upon 
Tl-fiUing accompagnied by a substitution of one Co atom 
per cell by Fe. 

Due to the negligible ionic character of the chemi- 
cal bonds in CoSba and TlFeCo3Sbi2, the static atomic 
charges of both compounds are rather small. In contrast, 
the dynamic Born effective charges are found to be very 
large and of the same amplitude as in ferroelectric oxides. 

The computation of the energy as a function of a co- 
operative displacement of Tl revealed that the center of 
the cage is the global minimum of the potential energy 
and not a saddle point as it has been suggested previ- 
ously. In the filled compound, the presence of Tl gives 
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FIG. 5: Phonon dispersion and projected DOS in CoSbs 
DOS show respectively Co/Fe, Tl and Sb vibrations. 



(a) and TlFeCo3Sbi2 (b). The fuU, dashed and dotted hnes in the 



rise to a single peak in the phonon DOS. Due to the hy- 
bridizations between the vibrations of Tl and Sb, this 
peak is significantly lower in frequency than the Tl bare 
frequency computed from the curvature of the potential 
energy well. However, these hybridizations are modest 
compared to those of the La-filled skutterudites where 
they give rise to a second La peak in the DOS. 

The Tl dominated mode at the F-point is polar and 
gives rise to a well defined band in the infrared reflec- 
tivity spectrum at low frequency and a strong increase 
of the static dielectric constant in the filled compound. 
On the one hand, this result shows that infrared spec- 
troscopy can act as a sensitive probe to study the lattice 
dynamics of the filling atom in skutterudites and related 
compounds such as clathrates. On the other hand, it 
suggests that the dielectric constant of materials charac- 



terized by an open structure containing empty voids can 
be tuned by filling the voids with foreign atoms. 
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FIG. 6: Tl projected DOS in TlFeCo3Sbi2 (full line), Tl 
projected DOS where we artificially decreased the mass of 
Tl to that of La (dashed line) and Tl projected DOS where 
we changed both the mass and the curvature of the potential 
energy well to those of La (dotted line). The vertical lines 
indicate the position of the corresponding bare frequencies. 
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FIG. 7: Theoretical and experimental heat capacity of CoSbs 
and Tl-filled CoSba divided by the temperature as a function 
of the square of the temperature. Experimental results are 
from Ref. H 
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